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Введение. Несоответствие структуры математической модели 
наблюдаемому процессу является одной из причин срыва 
сопровождения при функционировании управляемых систем. 
Это определяет актуальность задачи структурного синтеза 
следящих адаптивных систем, которая рассматривается в экс- 
тремальной постановке 

Материалы и методы. Для решения задачи синтеза в квази- 
детерминированном приближении используется методология 
объединенного принципа максимума. Исследование в этом 
случае заключается в поиске условия минимума целевого 
функционала, который является результатом свертки крите- 
рия качества и интеграла действия. Это позволяет построить 
модель фильтра состояния со структурной адаптацией как 
решение обратной задачи динамики. Практическая реализа- 
ция полученных с применением такого подхода результатов в 
цифровых следящих системах ставит задачу выбора схемы 
конечномерной аппроксимации полученной модели, что свя- 
зано с определением вида приближения синтезирующей 
функции и использованием общего решения однородного 
дифференциального уравнения. 

Результаты исследования. Рассмотрен вариант практической 
реализации цифровых следящих систем на базе конечномер- 
ной аппроксимации с новой структурой переходной матрицы 
и ковариационной матрицы вектора интенсивности возмуще- 
ний в дискретном времени. 

Обсуждения и заключение. Результаты математического мо- 
делирования показывают, что новый метод синтеза следящих 
систем в дискретном времени в условиях регулярных возму- 
щений обеспечивает повышение точности функционирования 
алгоритмов оценки в сравнении с традиционными. 


Ключевые слова: дискретное время, структурная адаптация, 
объединенный принцип максимума, синтез, фильтр Калмана. 
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Введение. Наиболее общий подход к синтезу следящих систем связан с решением стохастического диффе- 
ренциального уравнения Стратоновича. Его конечномерная аппроксимация приводит к сеточным фильтрам с боль- 
шим объемом вычислительных затрат. Поэтому практическое применение получил ставший эталоном для сравнения 
фильтр Калмана. Структурно подобные алгоритмы, например, расширенный фильтр Калмана, регрессионные филь- 
тры, сигма-точечные фильтры и т. д., применяют в связи с необходимостью оценки нелинейных систем в условиях 
воздействия негауссовских случайных процессов [1-4]. Эффективность таких решений не очень высока, поскольку 
они используют аппроксимацию, обоснованность применения которой не очевидна и сводится к эмпирическим пред- 
положениям. Применение метода расширения пространства состояний и построение модели наблюдаемого процесса 
на основе формирующего фильтра выглядит более обоснованно. Однако не всегда реальный нелинейный процесс 
можно аппроксимировать линейным фильтром с достаточной степенью точности. Это связано со структурными отли- 
чиями линейных и нелинейных моделей динамики, приводящими к возникновению ошибок в следящих системах, ко- 
торые можно трактовать как «структурные» и рассматривать задачу адаптации структуры модели с целью их компен- 
сации. При реализации фильтра Калмана эти ошибки частично компенсируются за счет возникновения корреляцион- 
ных связей, которые определяются решением уравнения для ковариационной матрицы ошибок фильтрации (оценки). 
Однако в установившемся режиме коэффициенты обратной связи подобных алгоритмов практически перестают зави- 
сеть от наблюдений и, соответственно, «структурные» ошибки полностью не компенсируются. На практике исполь- 
зуются менее точные, но более эффективные в смысле вычислительных затрат @а-В и а-В-у алгоритмы. Это позволяет 
сделать вывод о возрастании ошибок, связанных с неточностью применяемой в основе фильтра оценки состояния мо- 
дели [5, 6]. 

Применение квазидетерминистских методов в задачах синтеза следящих систем также приводит к калманов- 
ским алгоритмам. В частности, такой результат получают на основе методов теории регуляризации или процедуры 
инвариантного погружения [2,3,7]. Несоответствие положенной в основу процедуры синтеза математической модели 
наблюдаемому процессу в данном случае приводит к ошибкам экстраполяции и является одной из причин расхожде- 
ния фильтра и срыва вычислительной процедуры [2,4,6]. Попытки предотвращения такой ситуации связаны с исполь- 
зованием многомодельных алгоритмов или сложных моделей, характеризующихся большим числом подлежащих 
определению параметров [4]. Это приводит к возрастанию размерности задачи и препятствует реализации таких ре- 
шений в следящих системах в реальном масштабе времени из-за ограниченной вычислительной мощности электрон- 
ных вычислительных машин [1,9]. Таким образом, для повышения эффективности следящих систем необходимо ре- 
шить задачу адаптации структуры модели исследуемой системы [11,14]. 

Основой решения задачи структурной адаптации с успехом может служить методология объединенного 
принципа максимума [10-12]. Этот подход приводит к модели динамической системы, которая удовлетворяет прин- 
ципу Гамильтона-Остроградского [13]. Структура модели определяется из условия максимума функции обобщенной 
мощности с точностью до нелинейной синтезирующей функции, определяющей скорость диссипации и, соответ- 
ственно, степень структурной адаптации. Практическая реализация полученных с применением такого подхода ре- 
зультатов в цифровых следящих системах ставит задачу выбора схемы конечномерной аппроксимации полученной 
модели [9], что связано с определением вида приближения синтезирующей функции и использованием общего реше- 
ния однородного дифференциального уравнения. 

Целью данной работы являлась разработка метода синтеза следящих систем в дискретном времени с адапта- 
цией структуры модели на основе методологии объединенного принципа максимума. 

Постановка задачи. Пусть динамика системы в дискретном времени задана с точностью до структур пере- 


ходной матрицы Ф(К) = А””” и вектора интенсивности возмущений С (К) = А’ линейным разностным векторным 
уравнением [2,5,8] 

х(к+)=Ф()х()+6(Ющ®. (1) 
где х (К) = А’ вектор состояния, п(А ) —возбуждающая функция, К — момент времени. 

Уравнение наблюдений имеет вид 
У(к) = Нх(к)+Е(К), (2) 
где НЕА” — матрица проекции пространства состояний на пространство наблюдений, © (^ ) — белый гауссовский 
шум. 
В пространстве наблюдений выбран целевой функционал 


Ех 
= верни], (3) 


где К ' — параметр, который характеризует интенсивность помех в канале наблюдений, знак ^ означает оценку. 
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Требуется из условия минимума целевого функционала (3) найти оценки вектора состояния Х (К) = А’ для 
момента времени К ‚ если он увеличивается [3, 11]. 
Синтез следящей системы на основе объединенного принципа максимума. Для решения поставленной экстре- 


мальной задачи (1)-(3) необходимо определить структуру матрицы переходов состояний Ф(К) = А””” и вектора ин- 


тенсивности возмущений С(К) = А’ исходя из требований к синтезируемой следящей системе. 


Математическое описание динамики системы основано на физических закономерностях и допущениях и по- 
этому является, как правило, приближенным, что приводит к неадекватности модели, которая трактуется как струк- 


турная неопределенность Ф(К) = А”””. Например, так происходит когда в основу модели заложены лишь кинематиче- 


ские связи между параметрами движения [5,7]. 

Один из способов разрешения проблемы структурной неопределенности заключается в построении моделей с 
использованием утверждений, которые обеспечивают исчерпывающее описание динамики систем. За основу могут 
быть выбраны вариационные принципы [13]. В работе [11] с использованием принципа Гамильтона-Остроградского 
получена модель состояния в форме дифференциального уравнения второго порядка 














. в: 
= =е———+6 У(Ё)-Н ти й 
Ао) -н(я] ы 
где &= — параметр адаптации, определяющий долю регулярной составляющей приращения вектора состояния, 
1 ) . — 
с =—— —_ параметр, определяющий вклад случайного процесса в наблюдаемое изменение состояния, ОЕ Ср — 
^.Ке 
обобщенная сила, Со — замкнутое ограниченное множество. 
Модель (4) получена на основе условия [10] 
— у —1 ^ А . 
Фф = тах [АО (хх) — № [у(2) -Н (5,1) ]| х, ^ = СОИУЕ, (5) 
где Л. — неопределенный множитель Лагранжа, откуда с точностью до синтезирующей функции н(%%) [12] 
^ —1 л^ ^ А —1 ^ 
О=)^, + [(6)-н(%1] |. (6) 


Для получения конечномерной аппроксимации (4) в форме векторного разностного уравнения (1) запишем 
его как векторное дифференциальное уравнение первого порядка. Процедура расширения пространства состояний 
приводит к следующему результату 


0 
х= Л(х)+Ст(!), хи) =х , (7) 
х 
м А А - ы 0 1 0 
где х= 5х =Х, Х =Х, [ (х) = [х, |2 = А(х)х — вектор-функция системы для которой, А = ов а = = | 
х› -=——— 
№ 
. Е ре 
(2) — центрированный относительно Нх случайный процесс с интенсивностью К, ‚.В=-5 т 
т 


Рассматривается широко распространенный случай [5], когда выполняется с достаточной степенью точности 
условие: В = соп5Е . Тогда в дискретном времени для решения уравнения (7) справедливо [5] 


хе-+ АТ) = е АТ х(1) + р вт (г) ат, (8) 
[А 
откуда 
(А-ПАТ 
ФК) =е^^", СКи= |е^ АГ 9а пе, (9) 
КАТ 
где АТ — фиксированный интервал дискретизации. 


Вычисление матричной экспоненты приводит к матрице перехода состояний, которая определяется двумя па- 
раметрами 
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ВАТ _| 
Е 
Ф(АТ,В) =е^^Т = В , (10) 
0 с№Г 
а вектор интенсивности возмущений сверткой 
В 
(К+ПАТ е’ —1 (Е+ПАТ 
[10 п! (т) 
с(®®)= | В | кожи | | и о (11) 
т |0 ев о вт 12" 


В _ 
где пс, п) 6 


Решение полностью формализованной экстремальной задачи (1)—(3) с использованием дискретного варианта 
метода инвариантного погружения приводит к алгоритму калмановского типа [2-—4, 8] 


(к +1) = Ф(лт,в) (ЕК) 

(к) = (ик -1+к(ик-П(®-ня(-1] 

кк -1)= Рак Ин” [НРе-Инт +82". (12) 
Р(к—1)= Ф(лт,В)Р(к — А — 7 (Ат, В)+ У(К) 

Р(ик)= Р(ик- 1) (кант [НР +2 |" НРК -1) 


где я(Е+ИА) — оценка х(К-+1) по всем предыдущим измерениям, включая последнее  -е значение (предсказа- 
ние на шаг), #(ЕК) — оценка х() но всем предыдущим измерениям (фильтрация), К («К _ 1) — коэффициент об- 


ратной связи, Р( |) ‚Р (^ [А — 1) — соответственно ковариационные матрицы оценки и предсказания на шаг, 


СВ, [С -1 2-2еВАТ у СК, [7 1 т 
Е 
В 28” В? В В 2В В 
г. е2ВАТ _] 1-еРА7 Ги. е2ВАТ _1 
В | в В | | 2_ | 


—_ ковариационная матрица вектора возмущений. 





у — (13) 





р 


Исследование эффективности предлагаемого решения. Рассмотрим варианты построения следящих систем 
в ы 2 
с использованием предлагаемой модели и традиционно используемых вариантов кинематических моделей для хЕ^ 


ихЕА" [7] 








ж (А-ПАТ” 
где х› =х, х; =Х — соответственно скорость и ускорение. 
При 
м 
= ‚Н=П 0], (15) 
^2 
из (14) имеем [5] 
2 
Ф= - а = в 16 
ПВ т 
АТ 
откуда для вектора интенсивности возмущений 
2 
(А-ПАТ ЗАТ 
С(к)и(к)= | 2 |и(г)аг. (17) 


КАТ АТ 


Его использование позволяет получить ковариационную матрицу возмущений 
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И = (18) 
КАТ“ КАТ 
2 2 
Пусть 
ый 
х=| м |Н=1 0 0, (19) 
. 
где ускорение описывается моделью формирующего фильтра в непрерывном времени [5] 
х =-ах +, (20) 
я — постоянная времени ускорения, у’ — возбуждающая функция. Тогда при выборе возбуждающей функции в 


виде белого шума уравнения движения записываются в форме векторного дифференциального уравнения первого 
порядка, откуда [5] 
1 


1 АТ — [-1+аТ + ехр|- АТ] 
Ф=0 1 Тр — ехр|- АТ] (21) 
0 
0 ехр|- одТ| 


При фиксированном & и достаточно высокой частоте поступления измерительной информации, когда «АТ 


мало, (21) вырождается в матрицу Ньютона 


2 
Е В 
2 

Ф=б 1 ЛГ, (22) 
оо 1 


которая может быть получена из (14). Соответственно вектор интенсивности возмущений определяется выражением [5] 


Оки) = 


р Ат -т 1+ а + ПАГ + ехр[ о(& + ПА — <] 
(Е+ПАТ о, 




















1 23 
= | 0 1 — | _ ехр|- о(А + АТ —- т) О мс) ас. (2) 
(9 
о 0 ехр|- о(^к + АТ = т) 1 
Из (23) получены элементы ковариационной матрицы возбуждающей функции [5] 
2аК 
У = я т | -е ^^" +2аАТ +297 АТ” - далТе““" | 
(9) 
2аК 
= [-2аАт +1-е “^! +2адТе“" + 2аАТ + оАт? | 
24 
20 
Ма = 5 | НЫ. ст | 
24 
2аК 
у, [е-от 3 е-2аАТ + Залт| (24) 
3 
20, 
20 
о 204 а. г7чАТ | 
24 
_ 29 | _ -2оат | 
ы 24 
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Оценка эффективности предлагаемого решения проведена на основе сравнения результатов функционирования трех 
моделей (13), (18) и (24) [5,7]. 

Рассмотрен переходной режим функционирования в условиях шумов наблюдения, когда входное возмущение 
задано ступенчатой функцией. На рис. 1 представлена одна из реализаций переходного процесса. Здесь цифрами обо- 
значено: | — переходной процесс новой следящей системы (12); 2 — переходной процесс фильтра Калмана с моде- 
лью для равномерного движения; 3 — переходной процесс фильтра Калмана с моделью ускорения Зингера (21); 4 — 
возмущение. Итоги 100 вычислительных экспериментов подтверждают, что переходной процесс (12) выгодно отлича- 
ется меньшей колебательностью и временем сходимости от фильтров с (18) и (22). 


ея 
1.5103 ‚3 


о? о 


0 СоБыЫЫЕ || 60 70 
К 


Рис. 1. Переходной процесс 


Е1о.1. Тгапз1епе ргосез5 


Анализ результатов математического моделирования позволяет сделать вывод о том, что преимущество ново- 
го решения (12) в переходном режиме достигается за счет большего коэффициента обратной связи. Это демонстриру- 


ется на рис. 2. Здесь сохранены ранее принятые обозначения. 





0 20 40 60 К 


Рис. 2. Коэффициент обратной связи 


Е1с.2. ЕееДФасК рагатаег 


Выводы. Разработан новый метод синтеза дискретных следящих систем на основе объединенного принципа 
максимума. Он обеспечивает структурную адаптацию к наблюдаемому изменению состояния за счет новой линейной 
динамической модели в дискретном времени. Это проявляется в высокой скорости сходимости получаемых на его 
основе оценок к действительным значениям по сравнению с оценками фильтров Калмана с кинематическими моде- 


лями изменения состояния [5]. 
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